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Abstract 

The single-particle nuclear potential is intrinsically nonlocal. In this paper, 
we consider nonlocalities which arise from the many-body and fermionic na- 
ture of the nucleus. We investigate the effects of nonlocality in the nuclear 
potential by developing the Green's function for nonlocal potentials. The 
formal Green's function integral is solved analytically in two different limits 
of the wavelength as compared to the scale of nonlocality. Both results are 
studied in a quasi- free limit. The results illuminate some of the basic effects 
of nonlocality in the nuclear medium. 
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I. GENERAL ASPECTS OF NONLOCAL POTENTIALS 



A nonlocal term can be introduced to the non-relativistic Sclirodinger equation as follows: 

- |^V2M/(x) + [y^x) - E] M/(x) = -j rfVV^^(x,x') M/(x') . (1.1) 

Generally, such a nonlocal potential arises when one attempts to write the equation of 
motion for a single particle interacting with a many-body system. One source of nonlocality 
is the coupling of the elastic channel to other degrees of freedom. The potential in the elastic 
channel is modified by the contribution from excitation out of the elastic channel at a point x, 
propagation in an intermediate state, and then de-excitation back to the elastic channel at a 
point x'. This type of nonlocality is called Feshbach nonlocality. In addition, for a fermionic 
system, one must also account for exchange effects, which require antisymmetrization of the 
wave function between the projectile and the appropriate components of the target. This 
type of nonlocality is called Pauli nonlocality. 

The form of V(x, x') will not be considered in full mathematical generality. Neverthe- 
less, this does not seem to be a significant restriction for the physical theory. We consider 
the nonlocality to be characterized by a finite range h. In particular, we consider nonlocal 
potentials which limit to local potentials as the range 6 — 0. This form allows sufficient 
richness in the theory, but allows us to stay in contact with the familiar case of local quan- 
tum mechanics. This requirement excludes, for example, nonlocal potentials of the form 
Vi(x)V2(x'), which have no classical analog as a limiting case in the sense above ||l|. 

In the development below, we consider a nonlocal potential operator defined as 

Vnl = j I rf^xc^V ViVL(x,x') |x)(x'| . (1.2) 
For calculating the Green's function, we need the Fourier transform of V/vl(x, x'), i.e., 

V^l(p,p') = (p|VVl|p')- (1-3) 

Here and below a tilde is used to indicate the Fourier transform. Calculation of these 
matrix elements reveals some other general properties of the nonlocal potential. To begin 
with, consider the following form: 

Vnl (x, x') = 1^1 (x, x') 1^2 (x - x') . (1.4) 

We require that the second term limit to a delta function as the range of nonlocality 6 — > 0. 
That insures that the potential has a definite local limit. The first term is allowed for now 
to have an arbitrary dependence on x and x'. Then 

VW.(P.P')=(^//<i'xdV 

exp(-^) Vi(x,x')V2(x-x'). (1.5) 

We change variables to 



X exp 



2 



Y = ^ (1.6) 

and 

y = x-x' (1.7) 
so that the integrals can be partially separated: 

1 f j3 ^ ( ^(p + pO ■y\ 



VVl(p,p') = / d'yV2{y)ew 



X / d3Y\4(Y,y)exp f- '^^'f^"^ ) . (1.8) 



In principle, the second integrand depends on y. Since we assume that the potential is 
almost local, V2(y) is very sharply peaked in y. Then the y-dependence of Vi(Y,y) will 
not be seen at leading order in b. In that case, Vi{Y,y) ^ Vi(Y) alone. This is true to the 
extent that the range b of the nonlocality is small on the scale in which Vi is changing. For 
applications in nuclear physics, where the width b is of order 1 fm or less, this should be a 
very reasonable approximation. We have shown that the assumption of a small nonlocality 
range b motivates the following factorized form 

Vml{^, x') = Va (x - x') , (1.9) 

where Va (called Vi above) is the "average" local value of the nonlocal potential and 
(called V2 above) governs the "width" of the nonlocality. For a sufficiently small range of 
nonlocality, the particular functional form of the width term should not be too important. 
What is important is that some nonlocality is allowed, with a characteristic scale b. 

If this factorized form is assumed, then the integrals in Eq. (|1.8|) separate, and then one 
finds for the Fourier transform 

VNLip, P') = Ua (P - P') , (1-10) 

where, in the p-representation, Ua is the average potential term, and 11^ is the width term. 
These are given by: 

t/a = (X - X') (1.11) 

and 

f/. (P - P') = K . (1.12) 

This illuminates an interesting property of the factorized form. Under the Fourier transform, 
the "average" term of the potential in the x-representation becomes the "width" term of 



3 



the potential in the p-representation, and vice versa. In particular, there are the following 
correspondences (neglecting factors): 



X + x' < > P — p' 

X — x' < > p + p' . 



(1.13) 
(1.14) 



These relations depend only on the particular factorization of Vnl, and not upon the details 
of Va and Ko. 

Let us consider a few special cases. Consider first the width term Vw For a small 
range of nonlocality b, is a. narrow peak (say Gaussian) in the x-representation. Then 
in the p-representation, the average term Ua is a broad Gaussian. Thus nonlocality in 
the x-representation leads to a modified dispersion relation in the p-representation. The 
narrower the nonlocality, the flatter the dispersion relation. As 6 — 0, the free-particle 
dispersion relation is recovered. Now consider the average term Va. If it is constant, then 
its Fourier transform is proportional to a delta function, leading to a local potential in 
the p-representation. To the extent that it is non-constant, it leads to nonlocality in the 
p-representation (via 11^ = Va). 

Typically, a Gaussian form is taken for Vw, with V^ = K«(|x — x'|). Other forms could 
be chosen, but this seems adequate to describe the nuclear scattering data well 0]. Besides 
having an evident local limit, this particular factorized form meets some other requirements: 
Vnl should be symmetric under the exchange of x and x'; and when Va = const., as in infinite 
nuclear matter, V^l should be translation invariant. In this paper, we use the Gaussian form 
for Vw. In Fig. 1, a schematic matrix representation of the above discussion is given. 



As we have asserted in the previous section, nonlocality is a general characteristic of a 
single-particle description of a many-body system. It can be present, for example, in descrip- 
tions of systems in atomic physics [Q, condensed- matter physics 0], and quantum optics [^. 
Here, we are particularly interested in the specific nonlocal potentials which can appear in 
nuclear reactions. Below, we consider the nonlocality due to channel coupling (Feshbach 
nonlocality) and that due the fermionic character of the system (Pauli nonlocality). 

It is well known that a given nonlocal potential VnlI'^,^') can have a corresponding 
energy-dependent local equivalent potential Vl£;(x, i?) 0. Such a local equivalent potential 
gives the same cross section as the nonlocal potential. If the local equivalent potential also 
gives the same wave function, then it is said to be trivially equivalent. 

It was shown in a recent article [0 that most of the observed energy dependence in 
the inner region of the phenomenological local equivalent potential arises from the Pauli 
nonlocality. The Feshbach nonlocality would then be related to the energy dependence at 
the surface of that potential. Exchange effects are important at relatively short distances 
(inner region of the nuclear potential), where the density is high. Channel coupling effects 
are important at relatively large distances (nuclear surface), where the coupling form factors 
are peaked. As an example of this kind of long range effect we have, for example, the nuclear 
structure effects in Coulomb excitation. The relative contributions of the Feshbach and Pauli 
nonlocalities can also be studied using the technique of Ref . |Q] . 



II. NONLOCALITY IN NUCLEAR REACTIONS 
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We cannot give a rigorous proof that the factorized form for VArL(x, x') in Eq. ( |1.9| ) is 
completely general. However, we will motivate this form for both the Feshbach and Pauli 
nonlocalities below. 



A. Feshbach Nonlocahty 



Recent theoretical and experimental developments in the study of fusion reactions below 
the Coulomb barrier led to a renewed interest in the coupled-channels formalism |^. For a 
long time it has been known that coupling to nonelastic channels leads to a nonlocality in 
the elastic channel . The factorized nonlocal potential of Eq. ( p..9|) is sometimes adopted 
in channel coupling calculations, and seems to provide a good description of the energy 
dependence of the nuclear potential. The generation of a nonlocality in the elastic channel 
due to channel couplings has been described in detail in Ref. Here we consider simplified 
example, using only two channels. The coupled Schrodinger equations in this case are given 
by: 







2m 




+ Vi 




E 



*o(x) 
^i(x) 



(2.1) 



where \E'o(x) and ^l'i(x) are the wave functions of the ground and first excited states, re- 
spectively. Vl(x) is a local potential, -F(x) is the coupling potential and e is the excitation 
energy to the first state. These differential equations can be solved by the Green's function 
method. We can decouple these equations by eliminating either channel. We first solve for 
vl>i(x): 



vi/i(x) = j^j rfV G(x,x') F(x') vI/o(x') 
with the Green's function ^(x, x') is given by 



(2.2) 



G'(x,x') = (x 



E — /2m — Vl(x) — e ±iT] 



(2.3) 



Then the equation for \E'o(x) can be written as 



^o(x) = -J d^x' F(x) G(x,x') F(x') ^o(x') 



(2.4) 



Now that the \E'i channel has been eliminated, this is an equation for \E'o alone. By comparing 
Eq. ( p.4|) to the general form in Eq. ( |1 . 1|) , we find that the channel coupling generates a 
nonlocal potential 



r^L(x,x') = F(x) G(x,x') F(x') 



(2.5) 



The nonlocality in (^(x, x') arises because the kinetic energy is not diagonal in the x, x' 
representation. This and the elimination of all but the elastic channel leads to the nonlocality 
in the elastic channel. 
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As a simple example, we can assume that F{x.) has an approximately exponential shape 
and that Vi(x) has a slow variation. Then we can write Eq. ( p.5|) as (keeping in mind that 
in this situation we have (^(x, x') ~ (^(x — x') ||12|| ): 

yjvL(x,x') ^F(x + x') G(x-x'), (2.6) 
which has the same form as Eq. ( |1.9|) . 



B. Pauli Nonlocality 



The simple factorized form as given by Eq. (|1.9|) has a clear local limit when the nonlo- 
cality range b goes to zero. There are good reasons to consider such a simple form for the 
nonlocal potential V^ri;,(x, x'). For example, Bauhoff et al. |jl3| examined Pauli nonlocality 
in nucleon-nucleus interactions by using a folding approach. Adopting a nucleon-nucleon in- 
teraction written as the sum of a direct term and an exchange term, these authors obtained 
a nonlocal nucleon-nucleus potential from the microscopic description. Aoki and Horiuchi 
using the Resonating Group Method, verified that the microscopic interaction can be 



decomposed as a sum of direct and exchange parts. A comparison of the results for the 
trivially and local equivalent potentials and those potentials from the Perey-Buck analyses 
was made in Ref. |15|. It was shown that a large portion of the nonlocality present in the 



Perey-Buck potential is due to exchange effects. The microscopic approach of Ref. was 
found to produce very similar results to those using the Perey-Buck factorized form given 
above in Eq. (O). 



We conclude that the effects of the Pauli nonlocality in nucleon-nucleus interactions can 
be successfully described by a nonlocal potential given by Eq. ( |1.9D . The same conclusion 
can be reached in the case of nucleus-nucleus interactions |[7|JT6[]. 



C. Nonlocality and the Nuclear Dispersion Relation 



Here we note a consequence of allowing a nonlocal potential, as in Eq. ( |1 . 1| ) to describe 
the nuclear interaction For example, it is well known that the interaction described 

by Eq. ( |1 . 1| ) is dependent on the energy of the interacting system. In order to highlight the 
effects of nonlocality, we consider as an example an infinite nuclear medium (average term 
Va constant) and a Gaussian form for \4)(x — x'). Further, we neglect any isolated local 
potential Vi(x). We can then rewrite Eq. (|L9|) as 



V^7Vl(x,xO 



^3/25; 



exp 



52 



(2.7) 



Here h is the parameter describing the range of the nonlocal effects and Vq is the depth of 
the nuclear potential. In this way, we may rewrite Eq. ( p. . 1|) as 



■ ^2 r.' 

—V' + E 
2m 



IT 



3/2^3 



(i^x' exp 



62 



*(x')- 



(2.8) 



We first change the integration variable to s, defined by 6s 
exp(6s ■ V) ^(x) = \E'(x + 6s), we can write 



x' — X. Then using the relation 



6 



2m 



V' + E 



(2.9) 



We replace ^'(x) with its definition via the inverse Fourier transform: 

3/2 



(2.10) 



From this we obtain the momentum-space representation of Eq. ( |2.9[ ): 

.2 \ / p252\ 



|- + EU(p) = yoexp 



Ah' 



<|.(p). 



(2.11) 



Defining p = and k 
matter is given by 



|k|, the dispersion relation for a nucleon moving in infinite nuclear 



2l2 



2m 



E — Vq exp 



k'b 



(2.12) 



Equation ( p. 121) shows the dependence of the nuclear potential on momentum. Note that 
the strength of the potential decreases quickly with increasing momentum. For a local 
interaction (6 = 0), we recover the standard dispersion relation: 



2l2 



2m 



E-Vn. 



(2.13) 



Since Eq. ( |2.12| ) is a transcendental equation, it cannot solved for general momenta. For 
small momenta, the exponential can be expanded in a Taylor series. Below, only the two 
lowest orders of this expansion are kept. Thus we require that kb <^ 1, which means the 
nucleon wavelength is much larger than the range of the nonlocality. That is reasonable 
for the small range of nonlocality that we assume. Then the solution to Eq. ( |2.12D can be 
rewritten 



2m* 

where the effective mass m* is given by 



m 



{E - Vo) 



m 



1 - {mh'^Vo/2h') 



(2.14) 



(2.15) 



Note that we have changed k to k above to emphasize the fact that m* is used in the 
definition of the latter. Since Vq is negative, m* < m. The nonlocal interaction between the 
nucleon and the nuclear matter changes the dispersion relation, lowering the effective mass. 
In terms of the dispersion relation, a nucleon with a mass m and an nonlocal interaction is 
approximately equivalent to a nucleon with a mass m* and a local interaction. 
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III. DEVELOPMENT OF THE NONLOCAL GREEN'S FUNCTION 



A. General Development 



In this section, we study the Green's function for a nonlocal potential and present our 
main results. This allows a unifying approach to various limiting will be shown 

below. Studies of the general properties of nonlocal potentials and some limits thereof have 
been made in other formalisms |]T^,PD[ . These analytic studies are important in part because 
of the difficult nature of solving these systems numerically (but note a new approach in 
Ref. 0]). The Green's function may also be useful in calculating the transmission coefficient 
in heavy- ion collisions, where a nonlocal potential modifies the usual result ||22|| . 

We begin the development by defining the Green's operator for a nonlocal potential as 



G 



1 



NL 



E - Hnl ± iv 



(3.1) 



where 



H 



NL 



2m 



NL 



and Vnl is the nonlocal potential operator given by Eq. 
will be considered in a later section. 

The Green's function in the x-representation is given by 



(3.2) 

An additional local potential 



GArL(E;x/,x 
This can be rewritten as 

G'ArL(£';X/,X, 



2m ^ ^ 



Gnl I Xj) 



2m 



1 



E - Hnl ± iv 



(3.3) 



1' 
2m 



— / / d'p d^p' 



X (x/ I p) ( P 



1 



E - Hnl ± ir] 



P' (P' I X.) 



with 



P) 



exp(?p -^f/K) 



(3.4) 



(3.5) 



(27r^)3/2 ■ 

To proceed further with Eq. (|3.4| ), we adopt a factorized form for the nonlocal potential Vnl 
as in Eq. ( iLQl) . 

In order to understand the effects of nonlocality on the Green's function, we consider a 
simple nonlocal potential. We take the local or average term in this potential to be constant. 
For now, the form of the nonlocal or width term will be left unspecified. Thus 



VArL(x,x') = Va 



'x + x'' 



\/oK,(x-x') 



(3.6) 
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where Vq is a negative constant, and Vw is peaked at x = x'. As discussed above, taking 
Va constant forces f/^ = to be a delta function, so that the potential is local in the 
p-representation. That allows one of the integrals in each of Eq. ( p. .31 ) and Eq. ( p.4| ) to be 
collapsed. Then Eq. (|1.3|) can be written as 



(3.7) 
(3.8) 

can be written 
(3.9) 



(P I Vnl I p') = S'{p - p') Vo J d'y K,(y) exp [-z y ■ p/2h] 
where y = x — x'. Then 

1 A ^ '^^(p-pQ 

E-HNL±iri I E- pV2m - Vq K,(p) ± tv ' 



P 



where a tilde indicates the Fourier transform. Using these results, Eq. 



h r 

GNLiE;Xf,Xi) = — J (fp 



1 



exp[2p ■ (x/ - Xi)//i] 



{2Trh)^ E - p2/2m - Vq K,(p) ± iv 
For definiteness, we now choose for the nonlocal or width term 

= ^ exp (-^) . 



(3.10) 



As noted above, this form has a simple local limit when 6 — 0. This implies for the Fourier 
transform 



Vu,{p) = exp 



p^ 

Ah' 



(3.11) 



and thus for Eq. (|3^ 



1 



exp [ip ■ (x/ - Xi)//i] 



{2iTh)^ ^ - p2 /2m - Vo exp(-p262 /^n^) ± ' 
Defining new variables via E = h^k'/2m and p = hq (and q = |q|), this becomes 

I roo _ r+l 



(3.12) 



Gnl{.E] xj, Xj 



(277)2 7o 



X 



dqq J d{cos6) 
exp {iq\^f — Xj| cos( 



k^-q^ - (2mVo/n'^) exp{-qW/A) ± ir] 

Making the integration over cos^, we obtain 

1 i 



(3.13) 



GNL{E;:s.f,:s.i) 



^7r2 |x/ — Xj 



dq q 



(exp (ig|x/ - Xj|) - exp (-ig|xj - Xj|)) 
g2 + a exp(— /3g2) — (A;2 ± i^j^ 



(3.14) 



where 



9 



a 



2mVo 



(3.15) 



and 



P 



4 



(3.16) 



Equation (|3.14| ) is the starting point for the hmits of the subsequent sections. To make 
the integration over q in Eq. ( 3.14| ), we use the method of residues. This requires that we 
identify the poles of the integrand in the complex-g plane. Unfortunately, the equation 



q"^ + a exp - (fc^ ±ir]) = 



(3.17) 



is a transcendental equation, and the roots cannot be found in the general case. Below, we 
consider two limiting approximations involving (3 in which the roots of this equation can 



be found and the integral in Eq. ( 3.14 ) done analytically. In Fig. 2, the problem of finding 
these roots is indicated schematically. 



B. Digression on the Poles of the Integrand 

The integral in Eq. ( |3.14| ) is completely specified by the residues at the roots in q of 
Eq. ( p.l7|) . We show here how to develop expansions for these roots. In two cases of 



interest, these expansions can be truncated with minimal (and calculable) error. The errors 
on the integral results can be calculated using the known errors on the roots. Since the final 
results are based on controlled approximations, it is easy to determine where to truncate 
the expansions. 

Defining y = {q/kY, A = a/k'^, and B = pk'^, Eq. ( p.l7| ) can be written in a completely 
dimensionless form as 

y + Ae-^y -1 = 0. (3.18) 

For convenience, the ir] term has been dropped. From Fig. 2, one can see that there are 
exactly two roots on the real axis. We ignore any possible complex roots as they would 
not lead to pure propagating wave solutions in the Green's function. We exclude the case 
|y4| 3> 1; otherwise the value of A is unrestricted. 
Limit of E = /jfc^ < 1 : Rewrite Eq. (|3T7|) as 



y = l- Ae~^y 

y=(l-A) + A{By-]^B^y^ + . 

{l-AB)y={l-A) + A [-\bW + •.•)• (3.19) 

These are exact if all terms are kept. For i? -C 1, the latter term is small. The corrections 
can be solved for by perturbation, assuming the form 

(1 -AB)y = y^ + By, + B% . . . , (3.20) 
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and solving order by order. Then yo = 1 — A, and the first correction (the yi term) vanishes. 
Then Eq. ( 3.17| ) can be written as 

(1 - ap) q^-{e-a± ir]) + O [k^P^f) = . (3.21) 

This is a direct replacement for the denominator of Eq. ( p.l4| ) near the poles; note that we 
did not divide out the factor (1 — a/3). 

Limit of E = /jfc^ > 1 : Rewrite Eq. (|3T7|) as 

y = l- Ae-^y 

y=[l- Ae-^) - A (e-'^y - e"^) . (3.22) 

These are exact. For i? ^ 1, y ^ 1, so the latter term is small. The corrections can be 
solved for by iteration to yield a power series in . Note that B ^ 1 and ^ 1, so 
that Be~^ <^ 1. Then Eq. ( p.l7| ) can be written as 

g2 -(^e-a exp{-(3e) ±i7])+0 (k^ exp(-2/5A;2)) = . (3.23) 

This is a direct replacement for the denominator of Eq. ( |3.14| ) near the poles. 



C. Low-Momentum (or Effective Mass) Limit 

If (3k'^ <^ 1, then the denominator of of Eq. (|3.14|) can be rewritten as in Eq. (|3.21|) and 
so Eq. (|3.14|) becomes 

G7VL(£/;x/,Xi) = — / dqq 

OTT^ |Xj — Xj| m J-oo 

^ [exp (ig|xf - x^l) - exp (-ig|xj - x^|)] 



— (k^ ± irj) 



where 



is the effective mass and 



m 



m 1 — aj3 



(3.25) 



is a new wave number that is equal to the wave number of the nucleon in a constant local 
potential with the effective mass m*. By construction, these expressions agree with those 
found in Eqs. (|2.14|) and ( |2.15|) . Note that the corrections to Eq. ( |3.21|) have been dropped. 

The integral in Eq. ( p.24|) is done by the method of residues, with the contour closed in 
the upper-half or lower-half complex plane as necessary. We obtain the following expression 
for GiVL(-E; X/, Xj) in the limit -C 1: 

r (F-^ 1 exp(±m|xf -x,|) /m*\ 

G^,(E,x„x,) - ^——-^ [—^ . (3.27) 
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Nonlocality modifies the mass (making the term in parentheses different from unity) and 
modifies the wave number (which is also modified by the constant average term in the 
potential.) The local limit is immediately recovered as m* — > m and k ^ k. On the other 
hand, if additional orders in Eq. (|3.21| ) are retained, then the m*/m prefactor is unchanged 
and only n is modified. 

Note that this Green's function is an eigenstate of p^/h^ = — with eigenvalue k^. 
Using the definition of in Eq. i ^TI^ , fc2 = 2mE/h^ and the definitions in Eqs. QHTT^ ) and 



( P.16| ), this can be written as 

2m* 2m* 



E-Vo. (3.28) 



This check of the integration by direct differentiation returns the correct dispersion relation 
at the order of approximation used in this case. 



D. High-Momentum Limit 

If jSk'^ ^ 1, then the denominator of of Eq. (|3.14| ) can be rewritten as in Eq. (|3.23|) and 
so Eq. ( |3.14| ) becomes 



1 i r°° 
GNL{E;^f,^i) = —- / dqq 

OTT"^ |X/ — Xj| J-oo 



^ [exp (ig|xj - x,|) - exp (-ig|xy - x^|)] 



— (/c^ — a exp{—(3k'^) ± irj) 
This is again done by residues, and is found to be 



1 exp 


±i^Jk'^ — a exp (— 


/3P) 


X/ -x,|) 


An 


X/ -Xi 





Gml{E; x^, X,) = ^-^ ^ . (3.30) 



The small term exp(— /5/c^) damps the effect of the constant average potential. Other than 
the modifed wavenumber, the result in Eq. ( p.30| ) is exactly the usual free-particle Green's 
function. Note that k is the free-particle wave number, and is not modified by the constant 
average potential. If additional orders in Eq. (|3.23|) are retained, then the form is unchanged 
but the wavenumber is corrected by powers of exp(— /5/c^). 

Note that this Green's function is an eigenstate of p'^/h'^ = — with eigenvalue fc^ — 
a exp(— /5A;^). Using the definition k"^ = 2mE/h'^ and the definitions in Eqs. ( p.l5|) and 
( p.l6| ), this can be written as 

— = E-Vo exp(-mb'^E/2h'^) . (3.31) 
2m 

This check of the integration by direct differentiation returns the correct dispersion relation 
at the order of approximation used in this case. 
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E. Quasi-Free Limit 



Above, we simplified the integral by considering the limits oi B = f3k'^ <^ 1 and B = 
(3k'^ ^ 1, respectively. The value of A = a/k"^ was not specified, except for excluding the 
case \A\ ^ 1. The case |A| <^ 1 was allowed, and to consider its consequences it is sufficient 
to apply this limit to the previous results. It is not necessary to develop a new expansion 
for the roots or to do another integral. 

In the low-momentum case, we took B <^ 1 and the result was Eq. ( p.27|) . If \A\ <^ 1, 
then \A\B = \a\P -C 1, so that m*/m —>■ 1 and ^ k"^ — a. 



1 exp {±.i\/k'^ — a 


X/ -Xil) 


An 







GNL{E;Xf,Xi) ^ . . (3.32) 



This is of the same form as the free-particle Green's function, with the wave number modified 
by the constant average potential. This is a local limit, since the scale (3 has disappeared. 

In the high- momentum case, we took B ^ 1 and the result was Eq. ( p.30| ). In the case 
\A\ <C 1, the result is unchanged at the lowest nontrivial order: 



1 exp (±z^/A;2 — aexp (— Ix/- — x,;| I 
GME; x;, X.) - --^ ^ ^ ^ . (3.33) 

471- I 



This is also of the same form as the free-particle Green's function, again with the wave 
number modified by a term depending on the interaction. In this case, the nonlocality 
still appears as part of the interaction. Note that it would appear that these two Green's 
functions share a common limit if in the second we take jSk'^ <^ 1. Strictly speaking, this 
should not be allowed since the derivation assumed pk"^ 3> 1. 



F. Comparison to Data 



From Eqs. ( p.27| ), (|3.30|) , (|3.32| ), and (|3.33| ) we can see some of the basic effects of 



nonlocality in the nuclear medium. Nonlocality, in the low- momentum limit, changes the 
constant mass of the system to an effective mass, which depends on the range of nonlocality 
and the strength of the potential. This effect, in the limit \a\l3 -C 1, is negligible. In the high- 
momentum limit, nonlocality acts on the potential in a more complex way. However, this 
limit is satisfied in many cases of interest to nuclear physics ||2^. In this situation, the static 
strength of the potential changes, under the effects of the nonlocality, to a exponentially 
energy-dependent potential, e.g., the equivalent local potential. Chamon et al. |^ verified 



such a dependence of the equivalent local potential on the energy for elastic scattering of 
many different systems at high energies. 

We can define an energy-dependent local equivalent potential from Eq. ( |3.30| ) as 

y(E) = Voexp(-/?A;2) 
or, substituting the expressions for [3 and /c^, as 

V{E) = Voexpi-^E\ . (3.34) 
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As a test of our results, we compare Eq. (|3.34|) to the values of the inner region of the 
phenomenological potential extracted from data for the systems ^^C + ^^C and ^^C + 
^^O. In these comparisons, we adopted for the range of the nonlocality b = 0.14 fm and 
b = 0.11 fm for the first and second systems, respectively. These values of b are extracted 
from the data, and are in good agreement with those from a single-folding model. The 
values of Vq also come from the data. For more details about the data used, see Ref. and 
references therein. 

In the upper panel of Fig. 3, versus E for the system ^^C + ^^C is shown. The 

circles represent the values of the phenomenological potential at r = 4 fm extracted from 
the data and the solid line represents the result of Eq. ( |3.34| ). In the lower panel of Fig. 3 the 
same quantities for the system ^^C + ^^O is shown. For the ^^C + ^^C system, \a\P = 0.085; 
for the ^^C + ^^O system, \a\(3 = 0.075. In both cases there is an impressive concordance 
between the extracted and predicted values for the potential. 



G. Born Series to Include a Local Potential 



Referring back to Eq. ( |3.4|) , we note that in order to proceed we needed to collapse the 
double integral over p and p'. To do that, we took the interaction to be local in the p- 
representation, i.e., assumed a constant average potential in the x-representation. (In the 
usual development of the Green's function for a local interaction, any constant potential 
is defined away by shifting the energy.) One way to generalize the usual Green's function 
integral is to allow a constant potential with a nonlocality, as done here. This introduces a 
V^(p), and the remaining single integral can be done in some cases. 

To introduce an additional local potential, one can proceed as follows. We define a total 
Green's function by 

G = (3.35) 

E-Hnl-Vl± it, 

where 

p2 

Hnl = -^ + Vnl . (3.36) 
2m 



Here, Vatl is a nonlocal potential operator, as given by Eq. (|1.2|) , and is a local potential 
operator. Then 

G = Gnl + Gnl Vl G , (3.37) 

with 

Gnl = — -f- -— • (3.38) 

E - Hnl ± « V 

The Born series can be developed in the usual way: 

G = Gnl + Gnl Vl Gnl + Gnl Vl Gnl Vl Gnl + ■■■ ■ (3.39) 

This Born series for the operator G may or may not converge depending on the value of 
the energy of the system or on the strength of the local potential. The order by order 
contributions to the matrix element (xj | G | Xj) can then be developed once the potential 
is specified. 
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IV. CONCLUSIONS 



In this article we presented a Green's function approach to the study of nonlocal nuclear 
potentials. We showed that Green's functions can be used to investigate various limits in 
a unified way. In the low-momentum limit the effects of the nonlocal interaction can be 
accounted for with an effective mass and a wavenumber modified by the strength of the 
potential and range of the nonlocality. In the high-momentum limit the Green's function 
for the nonlocal potential also differs from the local one by a modified wavenumber. In 
this case, it is a function of the free wave number, i.e., the nonlocal potential is written in 
terms of an energy-dependent local equivalent potential. This behavior is consistent with 
the results recently obtained by Chamon, et al. p3| , who showed that as energy goes up the 
equivalent local potential for the nucleus-nucleus interactions become shallower. It would 
also be interesting to generalize the Green's function method to multidimensional systems 
where the effective mass becomes a tensor [^]. 

In the two main limits which we considered, we kept only the leading order results. 
However, in both cases we have shown how to easily generalize these results to any arbitrary 
order. 

The results in this paper are readily applicable to nuclear matter. We expect our for- 
malism to be useful to incorporate Pauli non-locality effects in the study of neutron stars. 
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Figure Captions 



FIG. 1. Schematic matrix representations of a nonlocal potential. The solid dots indicate 
nonzero matrix elements, and the vertical bars indicate the relative magnitude of those 
elements. In the matrix in the top left, only representative elements are indicated. In the 
top row, the potential has a constant local term multiplied by a Gaussian nonlocal term 
with a range b. On the left, the potential is represented in the x, x' basis; on the right, in 
the the p, p' basis. The two bases are related by a Fourier transform (FT). In the bottom 
row, we consider the local limit (6 = 0) of the top row. 

FIG. 2. The problem of finding the root of the dispersion relation (as a function of 
momentum q) is indicated schematically in the figure. The figure is symmetric about g = 0. 
The parabolic curve indicates the kinetic energy term /2m. The Gaussian curve indicates 
the nonlocal potential term Vq exp(— g^6^/4), with Vq < 0. The roots are at the points where 
the curves cross. When 6g <^ 1, the Gaussian is extremely broad and the roots are near 
q = ±^2'm{E + \VQ\)/h. When bq ^ 1, the Gaussian is very narrow and the roots are near 
q = ±^/2mE/h. 

FIG. 3. The energy dependence of the local equivalent potential |V^(-E')| at r = 4 fm in 
the range 10 < E/A < 200 MeV/nucleon. The circles are data, and the solid line is the 
result of Eq. ( |3.34|) . (See the text for details.) The upper panel is for ^^C + ^^C, and the 
lower panel is for ^^C + ^^O. 
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